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Abstract 

New nucleosynthesis bounds on the oscillation parameters of active neu- 
trinos mixed with a sterile one are derived for the non-resonant case. The 
controversy in the literature whether to use the annihilation rate or the total 
reaction rate for the estimates of sterile neutrino production is resolved in favor 
of the annihilation rate. In contrast to previous papers, the restrictions on os- 
cillations of electronic neutrinos are weaker than those of muonic and tauonic 
ones. 
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Neutrino oscillations in the early universe, especially if active neutrinos are mixed 



& \ with sterile ones, would have a noticeable impact on primordial nucleosynthesis Q and 

this permits to obtain interesting restrictions on the oscillation parameters. The re- 

<D ! _ 

suits very much depend upon a possible MSW resonance transition || in the primeval 

plasma. In the case of resonance (5m 2 < 0) the oscillations are much more efficient, 

neutrino spectrum can be strongly distorted, and the lepton asymmetry in the sector 

of active neutrinos can be enhanced by several orders of magnitude. However the 

calculations in this case are very complicated and controversial conclusions have been 

reached by different groups. For the discussion and the list of references see the recent 

papers f||, ||. In the non- resonant case the calculations are much simpler but there is 

also a disagreement between the papers || and all the subsequent works, see e.g. 0- 

|| . In the first papers it was assumed that the probability of production of sterile 

neutrinos, u s , is proportional to the rate of (inverse) annihilation of active neutrinos 

into light fermions in the plasma, jann, while in all other papers it was argued that 
^lso: ITEP, Bol. Cheremushkinskaya 25, Moscow 113259, Russia. 
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the rate of production of v s is proportional to a much larger total scattering rate, 
Itot = lei + "farm- Correspondingly the BBN (big bang nucleosynthesis) bounds on 
the oscillation parameters would be much more restrictive. 

It is shown here that these arguments are not correct and the probability of v s 
production is indeed proportional to the annihilation rate in agreement with ref. ||, 
though more precise calculations presented below result in stronger bounds than 
found in the papers || but still weaker than those obtained in the above quoted 
papers [§]-[§, where the probability of production was taken to be proportional to 
7tot. The simple argument showing that the effect of oscillations vanishes in the 
limit of jann = is that in this limit the total number density of active and sterile 
neutrinos is conserved, n a + n s = const, (see below eqs. (0,§)) and roughly speaking 
their total energy density remains the same as in the absence of oscillations. In 
fact the situation is somewhat more complicated because at an early hot stage the 
equilibrium with respect to annihilation was also reached, so the obtained expressions 
do not permit to take the limit of vanishing 7 ann . 

The difference between the approach of the present paper and all the other ones, 
where the BBN bounds have been derived, is that in those papers the breaking of 
coherence of the oscillations was described by the simplified anzats to the r.h.s. of 
the kinetic equations:: 



where curly brackets denote anti-commutator and g is the interaction matrix. In 
flavor basis it has the only nonzero entry in (a,a)-position: 



The coefficient 7 is determined by the rates of all neutrino reactions in the plasma; 
the calculations can be found in refs. ||, [7|, [Uj |3|. However this expression is not 
satisfactory for our purpose, in particular, because it does not conserve particle num- 
ber in the case of elastic scattering. To this end we need the exact equation for the 



i{g 2 ,P-P (eq) } 





(2) 
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density matrix of oscillating neutrinos derived in refs. |l|, PH . However the Fermi 
blocking factors will be neglected in what follows. The contribution to the r.h.s. of 
kinetic equations from elastic scattering of neutrinos on some other leptons, I, in the 
plasma is given by the anticommutators: 

{^jf) l = ~ {A2J2) (f l ^9 2 ,p(Pi)} - fi(p4){9 2 ,p(ps)}) (3) 

where A e \ is the amplitude of elastic scattering properly normalized to give a correct 
result for the diagonal matrix elements. It is assumed that the leptons I are not 
oscillating, otherwise, if the latter are oscillating neutrinos, the matrix structure of 
the result would be much more complicated but in the first approximation we can 
use the expression (|3|). The integration over momenta of all particles except for 1 are 
assumed, namely the following integration of the r.h.s. should be done: 

1 f d 3 p 2 d 3 p 2 d 3 p 2 4 4 

Wj (2n) 3 2E 2 (2n) 3 2E 2 (2n) 3 2E 2 {2rr)6 & + " *> " & 

The amplitude of elastic scattering with proper symmetrization factors can be taken 
from tables of ref. JTO . 



Neutrino annihilation is described similarly: 

= -i A lnJ A ) [{9,p(pi)gp(P2)} + {g,p(P2)gp(pi)} - fMfiipd] (5) 



It is convenient to introduce real and imaginary parts of the non-diagonal compo- 
nents of neutrino density matrix: 

Pas= P* sa = R+H (6) 

where a and s mean respectively "active" and "sterile". 

Now the kinetic equations describing evolution of density matrix of oscillating 
neutrinos can be written as: 

paaiPl) = -FI - A 2 el [p aa (pi)fi(p 2 ) - paa{Pz)fl{Pi)] 

- A lnn [Paa(Pl)Paa(P2) ~ flfai) ffoi)] , (7) 



p„(pi) = FI, (8) 
= WI-{A 2 j2)[R{ Pl )f l (p 2 )-R(p 3 )f l (p 4 )) 

-(^Ln/4) [poo(pi)^(P2) + Paa(p2)R(pi)} , (9) 

j(pi) = -WR - (F/2) ( Pss - Paa ) - (A 2 J2) [I( Pl )fi(p2) 

-I(p3)fl(pd\ - (4L/4) [Paa(Pl)/(p 2 ) + P«a(P2)/(Pl)] • (10) 

Here we use notations of ref. [[J, so that 

F = 5m 2 sin 29 /2E (11) 

and 

W = 5m 2 cos20/2E + C l (G 2 F T i E/a) (12) 
where a = 1/137, Gf = 1.166 • 10~ 5 GeV~ 2 , and the constant C\ depends upon the 



neutrino flavor, C e = 0.61 and C MiT = 0.17 fpL2|| . We neglected the term related to 
charge asymmetry of the plasma, it may be a good approximation in non-resonant 
case. 

In the cosmological FRW background the time derivative in the l.h.s. of kinetic 
equations goes into d/dt — > d/dt — Hpd/dp, where H is the Hubble parameter ex- 
pressed through the thermal energy density as 

3H 2 m 2 Pl = Tr 2 g* 4 

8vr 30 1 ' 

The factor g* = 10.75 is the number of relativistic species in the cosmic plasma. 
We will introduce the new variables: 

x = m /T and y = p/T , (14) 

where the dimensional normalization factor mo is taken equal to 1 MeV. It is a good 
approximation to assume that the temperature evolves as the inverse cosmic scale 
factor, T ~ l/a(t). In other words, T = —HT and thus the differential operator 
(dt — Hpdp) transforms into Hxd x . 



Since the oscillation rate is much faster than the reaction rate (in other words, 
the terms ~ F, W are typically larger than the terms related to reactions) the equa- 
tion ([[!]) can be solved as 

F 



R 



2W 



Pa 



(15) 



This permits to express the imaginary part I through p aa using eq. ([5]) and to sub- 
stitute the result into eq. (^). If we assume that p ss is small in comparison with p aa 
so it may be neglected, the eq. ([?[) becomes a closed equation for a single unknown 
function p aa . If elastic scattering of active neutrinos is sufficiently strong (its rate is 
approximately an order of magnitude larger than the rate of annihilation), then one 
may assume that the distribution of active neutrinos is close to kinetic equilibrium 
with an effective chemical potential. In other words, the anzats p aa = C(x) exp(-y) 
is a good approximation. In this case both sides of eq. (0) can be integrated over d 3 y 
so that the contribution of elastic scattering disappears and the following ordinary 
differential equation describing evolution of C{x) is obtained: 

C 6 + 16(^ + ^/3 



dC 
dx 



Ki 

x i 



(~<2 

C 2 -l + — (6/ 2 
288 



ii 



96 1 + 2g\ + 2g% 



6J 2 - If 



where the constant K% is given by 



K, 



8Gl(l + 2gl + 2gl) 



7T° 



(16) 



(17) 



and 



/•oo r h \ roc 



dyy 3 e 



with 



2.34- 10" 



and /3 Mj . 



5m 2 cos 26 
see eqs. (0,0). 

The coupling constants of neutrinos to electrons are 



0.65- 10" 8 
5m 2 cos 29 ' 



(19) 



g\ = (0.5 ± sin 2 9w) 2 and g\ = sin 4 9 



w 



(20) 



The signs "±" refer to v e and v^ r respectively. With sin 2 9w = 0.23 we obtain 
9l + 9r = 0.5858 for u e and g 2 L + g\ — 0.1258 for v^ T . Correspondingly K e = 0.17 
and K VjT = 0.098. We assumed that F/W < 1 and thus the term ~ (F/W) 2 dC/dx 
was neglected. It is a good approximation even for not very weak mixing. 
Eq. ( |16D can be solved analytically if |5| = |1 — C| <C 1: 
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dx\Di(xi) exp 



2K ( 1 



(21) 



where 



D(x) 



Kx~ 



(6/ 2 + If) + k (6/ 2 - 1\ 



(22) 



with a = 3.47- 10 , b e = 4.37-10 , and b^ T = 5.55-10 . This expression determines 
the freezing temperature of annihilation of active neutrinos, TJ = (3/27^) 1//3 MeV. 
The increase of the total number density of oscillating neutrinos 

d 3 y 



i 727p (Paa + Pss) 



(27T) 

can be found from the sum of equations (|7j) and (||) and is given by 

I = / dx Di(x) 



n, 



e i J i 



(23) 



(24) 



All integrals can be taken analytically, first over dx and then over dy. Finally we 
obtain: 

2 



'An\ ttKi /sin 20 N 



[a, + 6, + (5.14/6)(o,-6,)] 



(25) 



This number should be smaller than the upper limit for extra neutrino species, AN V , 
permitted by BBN. In particular, for electronic neutrinos we obtain (for small mixing 
angles) : 



5m 2 sin 4 2B\ Ve < 5 • 10~ 4 AN 2 , 



(26) 



while for z/„ r the result is surprisingly stronger: 



5m 2 sin 4 201^ < 3.3 • 10" 4 AiV 2 



(27) 



The result makes sense only if AN < 1. Otherwise, even if the excitation of sterile 
neutrinos is complete, i.e. they reached equilibrium abundance, their maximum con- 



tribution to the effective number of neutrinos would be unity. The bounds ([26],|27|) are 
approximately an order of magnitude weaker than those obtained in ref. JFj and an 
order of magnitude stronger than found in ref. ||. In the previous papers the bounds 
on oscillation parameters of v e was stronger than those for and v T and it was re- 
lated to a faster production of u e by charged currents. However, as we saw above, 
there is a competing contribution related to the refraction index of active neutrinos: 
it is larger for u e , see eqs. (0101) ■ Hence the oscillations of z^ jT are less suppressed 
at high temperatures and their sterile partners are more efficiently produced. 

We neglected here a distortion of neutrino spectrum by the oscillations. In con- 



trast to the limits (0,^3) that are effective even for very small mixings, because the 
smallness of mixing can be compensated by the efficient production of v s at the early 
stages, the distortion of spectrum of v e could only be developed sufficiently late, at 
T < 2 MeV and the effect could be essential for sufficiently large sin 9. One possible 
form of spectrum distortion is a generation of an effective chemical potential of the 
same sign for v and u, while kinetic equilibrium shape of the distribution survives. 
This effect dominates for sufficiently large 5m 2 . It was estimated in ref. |J. A de- 
viation of electronic neutrinos from kinetic equilibrium is essential for smaller mass 
difference. It was accurately calculated for 5m 2 < 10 _7 eV 2 in the papers Spec- 
tral distortion may have a very strong effect on primordial abundances, so the limits 
on the oscillation parameters could be sensible even if the data permits AN V > 1. 

Electronic neutrinos are in a good thermal contact with the rest of the cosmic 
plasma till T = 2 MeVQ. Above this temperature one may assume that the spectrum 
of v e is close to the equilibrium Fermi-Dirac form (in fact we will take Boltzmann 
approximation, f u = f( eq ) exp(— E/T)). Below 2 MeV u e may be considered as free 
and the impact of medium enters only through the refraction index. One can check 

2 The reaction rate is proportional to neutrino momentum p, see e.g. ref. || where this issue is 
discussed. Thus the equilibrium is maintained down to a smaller T for larger p and vice versa. 
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that for 5m 2 > 5 ■ 10~ 6 eV 2 the vacuum term in effective potential ([T2] ) dominates and 
neutrino oscillations are close to the vacuum ones. The number density of electronic 
neutrinos in this case is equal to: 

p„.,„ a = (c 4 + s 4 + 2c 2 s 2 cos 5m 2 t/2E) f {eq \E), (28) 

where s = sin26 l and c = cos 29. The last term in this expression quickly oscillates 
at nucleosynthesis time scale, so it can be neglected. This term would give a non- 
negligible contribution to neutrino spectrum for 5m 2 < 10~ 8 eV 2 , but in this case the 
matter effects should be taken into account. From eq. (p§[) follows that the effective 
chemical potential of u e is 

£ = /VT = ln(l- S 2 /2) (29) 

As is argued in ref. || such a distortion of equilibrium neutrino spectrum is equivalent 
to a renormalization of the Fermi coupling constant Gp — > Gp[l + exp(£)]/2. Since 
£ < the freezing of n/p-ratio would take place at higher T and more 4 He would be 
produced. This confines the mixing of u e with v s to the region: 

sin 2 29 < 0.32 AN„ (30) 

for all 5m 2 > 5 • 10~ 6 eV 2 . 

For smaller mass differences matter effects cannot be neglected and the effective 
mixing angle would depend on neutrino energy: 

• oa sin 2^ 

Sm 2&eff = l + 1.22E*T*(Gp 2 /a)\5m 2 \-i [d } 

This would modify the energy spectrum of v e and change the frozen value of neutron- 
to-proton ratio. Detailed calculations for \5m 2 \ < 10~ 7 eV 2 can be found in the 
HT3I. 



papers 



The observational limit on AN V was analyzed in a recent review [14| with the 



conclusion that AN V < 0.2. Thus a large mixing of with a sterile partner proposed 



for explanation of the atmospheric ^-deficit [15| is excluded in non-resonant case. 



However, as is argued in ref. Hlq] , one should be cautious in the data analysis and 
probably the safer limit is AN U < 1. If this is true an arbitrary strong mixing of active 
neutrinos and v T with one sterile companion would be permitted by BBN. In the 
case of f e -oscillations the distortion of v e spectrum could be essential and in this case 
an interesting bound can be observed even if AN U = 1 is permitted. However if one 
admits existence of one sterile neutrino, it is natural to assume that there are three of 
them, as e.g. in the case of Dirac- Major ana mass mixing [lj]. In this case the mixing 
angles should be rather strongly bound by BBN, but detailed analysis with correct 
description of decoherence effects is necessary. 
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